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ABSTRACT 

We present the results of the timing analysis of five Rossi X-ray Timing Explorer observations 
of the Black Hole Candidate GRS 1915+105 between 1996 September and 1997 December. 
The aim was to investigate the possible presence of a type-B quasi-periodic oscillation (QPO). 
Since in other systems this QPO is found to appear during spectral transitions from Hard to 
Soft states, we analyzed observations characterized by a fast and strong variability, in order 
to have a large number of transitions. In GRS 1915+105, transitions occur on very short time 
scales (~ sec): to single them out we averaged Power Density Spectra following the regular 
path covered by the source on a 3D Hardness-Hardness-Intensity Diagram. We identified both 
the type-C and the type-B quasi-periodic oscillations (QPOs): this is the first detection of 
a type-B QPO in GRS 1915+105. As the spectral transitions have been associated to the 
emission and coUimation of relativistic radio-jets, their presence in the prototypical galactic 
jet source strengthens this connection. 

Key words: accretion, accretion discs - X-rays: binaries. 
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1 INTRODUCTION 

Systematic variations in the energy spectra and intensity of tran- 
sient Black-Hole Candidates (BHC) have been recently identified 
in terms of the pattern described in an X-ray Hardness-Intensity di- 
agram (HID, see Homan et al. 2001, Homan et al. 2005b, Belloni 
et al. 2005). Four main bright states (in addition to the quiescent 
state) have been found to correspond to different branches/areas 
of a square-like HID pattern. In this framework much importance 
is given to the intermediate states (called Hard Intermediate State, 
HIMS, and Soft Intermediate State, SIMS) and to the transitions 
between them, identified from the behaviour in several bands of 
the electromagnetic spectrum (from radio to hard X-rays, see also 
Fender, Belloni & Gallo 2004 and Homan et al. 2005b) and from 
the timing properties of the X-ray light curve. 

Low-frequency Quasi-Periodic Oscillations (LFQPOs) with 
centroid frequency ranging from mHz to tens of Hz have been ob- 
served in the X-ray flux of many galactic BHCs since the '80s (see 
van der Klis 2006; McClintock & Remillard 2006 and references 
therein). Three main types of LFQPOs, dubbed Type- A, -B and -C 
respectively, originally identified in the light curve of XTE 11550- 
564 (Wijnands et al. 1999; Remillard et al. 2002), have been seen 
in several sources (see Casella et al. 2005 and references therein). 
We summarize their properties in Table[T] 
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In the context of the state classification outlined above, it is 
possible to ascribe the three LFQPOs to different spectral condi- 
tions (see Table [T] Homan et al. 2001, Homan & Belloni 2005, 
Belloni et al. 2005). 

The type-C QPO is associated to the (radio loud) HIMS and 
to the low/hard state. It is a common QPO seen in almost all 
BHCs with a variable centroid frequency correlated with the count 
rate, a high fractional variability and a high coherence (Q = 
zy/FWHM~10). 

The type-B QPO has been seen only in few systems, although 
it is being seen in a growing number of sources (see Casella et al. 
2005 and references therein). It is a transient QPO associated to 
spectral transitions from the (radio loud) HIMS to the (radio quiet) 
SIMS. Its features are a ~ fixed centroid frequency (around ~6 
Hz), lower fractional variability and Q than type-C. Some authors 
(Fender, Belloni & Gallo 2004; Casella et al. 2004) suggested that 
these spectral transitions are in turn associated to the emission and 
coUimation of transient superluminal relativistic jets visible in radio 
band. These jets are seen in a number of sources (GRS 1915+105, 
XTE J1550-564,GX 339-4, XTE J1859+226, GRO J 1655-40, etc.). 
However, not in all of these sources we could resolve the spectral 
transition to see the transient QPO. 

The spectral properties connected to type-A QPO are similar 
to those introduced for the type-B. This QPO has been seen in few 
systems (Casella et al. 2005). It is broader, weaker and less coherent 
than the type-B QPO. 
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Table 1. Summary of type-A, -B and -C LFQPOs properties (from Casella 
et al. 2005). 



Properties 


Type-C 


Type-B 


Type-A 


Frequency (Hz) 


~ 0.1 - 15 


~ 5 - 6 


~ 8 


Q(;^/FWHM) 


~ 7 - 12 


Z 6 


S 3 


Amplitude (% nns) 


3-16 


~ 2 - 4 


S 3 


Noise 


strong fiat-top 


weak red 


weak red 


Pliase lag" <SvQPO 


soft/hard'' 


hard 


soft 


Phase lag @2uqpq 


hard 


soft 




Phase lag @i^Qpo /2 


soft 


soft 





" With "hard lag" we mean that hard variability lags the soft one. 
^ Trend towards soft lags for increasing QPO frequencies 

GRS 1915+105 is a transient BHC discovered on August 15 
1992 with the WATCH instrument on board GRANAT (Castro- 
Tirado et al. 1992, 1994). It is the first galactic source observed 
to have apparently superluminal transient relativistic radio jets 
(Mirabel & Rodriguez 1994), commonly interpreted as ejection of 
ultra-relativistic plasma, with a speed close to the speed of light (up 
to ~ 98%). Its radio variability was discovered to correlate with the 
hard X-ray flux (Mirabel et al. 1994). Thanks to VLA-radio obser- 
vations of these jets, Rodriguez et al. (1995) estimated a distance 
of 12.5 kpc; more recent estimates attest a distance of 6. 5 ± 1.6 kpc 
(Kaiser et al. 2005). The mass of the compact object was estimated 
through IR spectroscopic studies (Greiner, Cuby & McCaughrean 
2001, Harlaftis & Greiner 2004) to be 14.0 ± 4Mq which unam- 
biguously makes GRS 1915+105 a BHC. 

The various and rich phenomenology of this source was clas- 
sified by Belloni et al. (2000): they analyzed 163 RXTE observa- 
tions, showing that the complex behaviour of GRS 1915+105 can 
be described in terms of spectral transitions between three basic 
states, A, B and C (not to be confused with the name of the LFQ- 
POs introduced before), that give rise to 12 variability classes. The 
non standard behaviour of GRS 1915+105 (it is a very bright tran- 
sient source continuing the same outburst started in 1992) was in- 
terpreted as that of a source that spends all its time in Intermediate 
States (both in its hard and soft flavors), never reaching the LS or 
the quiescence (see e.g. Fender & Belloni, 2004). 

GRS 1915+105 also shows strong time variability on time 
scales of fractions of second, revealing low- and high-frequency 
QPOs (see Morgan et al. 1997) whose properties (frequency and 
fractional variability) are tightly correlated with the spectral param- 
eters (Morgan et al. 1997; Muno et al. 1999; Markwardt et al. 1999; 
Rodriguez et al. 2002a, 2002b; Vignarca et al. 2003). In particular, 
all LFQPOs observed from this system can be classified as type-C 
QPOs. Although GRS 1915+105 makes a large number of fast state 
transitions, which have been positively associated to radio activity 
and jet ejections, no type-B QPO has been observed to date. 

In this paper we present the discovery with RXTE of the type- 
B QPO in the X-ray light curve of GRS 1915+105. The QPO was 
present during fast spectral transitions that we identify with the 
HIMS to SIMS transition observed in other BHCs. 



2 OBSERVATIONS AND DATA ANALYSIS 

We analyzed five RXTE/PCA observations GRS 1915+105 col- 
lected during AOl and A02, between 1996 September 22 and 1997 
December 17 (see Table |2}. We chose observations belonging to 
variability classes /i and P (according to the classification of Bel- 
loni et al. 2000, see Figure [T] for two examples of light curves). 
Restriction to AOl and A02 is due to three reasons: first, AOl 




Figure 1. 2-13 keV light curves for the two variability classes examined 
here. Bin size is 1 second. 



and A02 observations are those for which a thorough classifica- 
tion was made by Belloni et al. (2000), which removes the need for 
a complete classification of observations throughout the archive, 
a task beyond the scope of the present paper; second, AOl and 
A02 were in average characterized by a higher number of working 
PCUs compared to more recent epochs (this improves the statis- 
tics) and third, using data from different epochs, you get different 
hardness-intensity and color-color diagrams. We chose class ^ and 
P observations because they correspond to intervals when frequent 
fast transitions among the three spectral states are observed. 

The rationale behind this choice was to analyze a large num- 
ber of spectral transitions from hard (C) to soft states (A and B), 
when a type-B QPO would be expected. The transitions being fast, 
the QPO would not be detected for a single transition, but adding 
a large number of data corresponding to the same transition, the 
signal could become detectable. The total amount of usable expo- 
sure was 44648 seconds. For each observation we produced light 
curves in the PCA energy channels 0-35 (corresponding to the en- 
ergy range 2-13 keV) and two color curves defined as the ratio be- 
tween the counts in the 5-13 and 2-5 keV bands (HRl, absolute 
channels 14-35 and 0-13) and between the counts in the 13-38 and 
2-5 keV bands (HR2, absolute channels 36-103 and 0-13). The bin 
size was 2 seconds for all these curves. Given the high count rates 
we did not subtract the background. In addition, for each 2-second 
time interval we accumulated a power density spectrum (PDS) in 
the 2-5 keV, 5-13 keV and 13-38 keV energy ranges. All the PDS 
were extracted with a time resolution of 1/128 s (corresponding to 
a Nyquist frequency of 64 Hz). The PCA data modes used to pro- 
duce the PDS were Binned and Event Data. For an explanation of 
the different PCA modes see for example Jahoda et al. (1996). 

For each observation, we created a 3-dimensional hardness- 
hardness-intensity diagram (HHID), where the source follows 
a regular path. For each observation we averaged PDS on a 
variable number of points (~ 100-800, to have enough statis- 
tics), following the whole diagram. The selections were made 
using the free software for multi-variable management XGOBl 
ii.http://www.research.att.com/areas/stat/xgobi/). For each selection 
we rebinned the average PDS logarithmically, and the Poissonian 
noise, including the Very Large Events (VLE) contribution (Zhang 
1995, Zhang et al 1995), was subtracted. The PDS were normal- 
ized to square fractional nns (see Belloni & Hasinger 1990) and 
fitted with a combination of Lorentzians (see Nowak 2000; Belloni 
et al. 2002). In a few cases (all the PDS of observation #1 and one 
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Table 2. Log of the 5 RXTE/PCA observations analyzed in this work 



N° 


Obs. Id. 


Date 


Starting MJD 


Exp. (s) 


Classification 


1 


10408-01-35-00 


1996 Sep 09 


50348.271 


9448 




2 


20402-01-45-03 (#1, 2, 3) 


1997 Sep 09 


57700.250 


10038 


P 


3 


20402-01-53-00 


1997 Oct 31 


50752.013 


9656 


/3 


4 


20402-01-53-01 & 20402-01-53-02(#l) 


1997 Nov 04-05 


50756.412 


6322 




5 


20402-01-59-00 


1997 Dec 17 


50799.091 


9784 


/3 



Obs. 4 is composed of two orbits from two separate observations 
# indicates the number of the RXTE orbit within the observation 
Classification is from Belloni et al. (2000). 




0,05 0.10 0,15 0,20 0,25 0.05 O.IO 0,15 0,20 0,25 
hr2 hr2 



Figure 2. Two-dimensional projections of the hardness-hardness-intensity diagram (HHID) for observation #5. Left side panel: color-color diagram; right side 
panel: hardness-intensity diagram. Definitions for the colors can be found in [|2] the count rate is the 2-13 keV energy band. The an'ows follow the source's 
regular and repetitive movement. Numbers correspond to the regions discussed in the text (see Figure[3). 



PDS of observation #3) the Poissonian subtraction left flat residu- 
als at high frequencies, hence we decided to account for it adding 
an additive constant to the model. The fitting was carried out with 
the standard XSPEC vlL3 fitting package, by using a one-to-one 
energy-frequency conversion and a unit response. 

For each 2-s interval, we also produced a cross-spectrum be- 
tween the 2-4.5 keV (absolute channels 0-11) and 4.5-1 1 keV (ab- 
solute channels 12-29) light curves. We then calculated averaged 
cross-spectrum vectors for all the different selected regions of each 
observation, from which we derived phase lag spectra (for details 
on the phase lag analysis see e.g. Casella et al. 2004). We defined 
phase lags as positive when the hard X-ray variability /oZ/ow^ the 
soft one. In order to quantify the phase-lag information for each de- 
tected QPO component we integrated the lags in a range centered 
on the QPO centroid frequency with a width equal to the FWHM of 
the QPO peak itself (see Reig et al. 2000). Notice that the estimate 
of the QPO lags depends on the relative power of the continuum 
and of the QPO itself. When the continuum power contribution is 



significantly larger than the QPO power the lags of the continuum 
may dominate. 



3 RESULTS 

An example of HHID is shown in the two-dimensional projections 
in Fig. |2] (from observation #5, class f3). The regular and repeti- 
tive movement followed by the source is marked by the arrows. 
Starting from the bottom-right corner of both projections, the PDS 
properties were found to change smoothly thorough the diagram. 
We could identify five main average behaviours of the PDS corre- 
sponding to five regions in the diagram. In Figure [3] we show the 
five regions in the two-dimensional projections of the HHID of ob- 
servation #5, together with five corresponding examples of PDS. 
In Figure |4] we show the evolution of the frequencies of all sig- 
nificant peaks (with significance a ^ 3) for all the observations, 
moving through the HHID following the arrows in Figure |2l The 
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Figure 3. Left panels: Two-dimensional projections of tlie HHID (see Fig. |2j- The tick points mark tlie five main regions. Right panel: Examples of five 
con'esponding typical power density spectra. For each power spectrum we also plot the best-fit Lorentzian components corresponding to the significant 
QPOs/bumps (ct 3). Dashed lines: type-B QPOs, continuous lines: type-C QPOs, dotted lines: harmonics. In the PDS in region 4 the Lorentzian fitting a 
bump identified as the harmonic of the type-B QPO (v = 13.84 Hz) has not been reported since it has been detected with a significance cr < 3. 



X axis has been chosen in such a way to have an equal spacing 
between different selections in the same region, since the same re- 
gion in different observations has been divided in a different num- 
ber of selections. This point can be clarified by looking tables IaTI 
to |A5| in appendix |A] where we reported all the significant peaks 
different from a continuum fitted in all the selections of the five 
observations. For each selection (in figure |4ll we reported the peak 
frequency obtained from the energy range where the feature was 
more significant. In four cases we found that the same feature was 
detected at more than 3 sigma as a QPO (i.e. with a quality fac- 
tor Q ^ 2) in one energy band and, with a higher significance, as 
a bump (Q < 2) in another. In these cases we decided to report 
the frequency obtained from the QPO. We considered two peaks as 
harmonics of each other if (i) their centroid frequencies, fitted in- 



dependently, are consistent with being harmonically related or (ii) 
fixing their centroid frequency to an harmonic ratio we could ob- 
tain a very good fit. We will no longer consider harmonic peaks in 
our discussion. 

In Figure |3] we can easily identify the presence of two peaks 
evolving in frequency thorough the diagram from region 1 to region 
4 (see Table[3]for a summary of their properties): 

(i) the first peak is detected (superimposed to a broad-band 
noise) since region 1. It increases in frequency (from ~2 Hz to 
~15 Hz) and quality factor Q (from 0.64 to 7.39) and decreases in 
fractional rms amplitude (from 22.86 to 1.16) while moving thor- 
ough the diagram from region 7 to 2 until it disappears at high count 
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Figure 4. Evolution of QPOs/bumps (with significance a ^ S) centroid 
frequencies as a function of tlie position in the HHID in Figure|3] for all the 
analyzed observations. Adopted criteria are explained in the text. 



Table 3. QPO types detected in the five regions for all considered observa- 
tions. Region 1 is not present in observations # 2 and 3. We indicated with 
"Bump-B" peaks with Q < 2. 



Obs. # 1 


Rej 

2 


;ions in 
3 


figure [5] 
4 


5 


1 


Bump-B, C 


C, B 


B 


No peaks 


2 C 


C 


C, B 


No peaks 


No peaks 


3 C 


Bump-B, C 


C, B 


No peaks 


No peaks 


4 


C 


C, B 


B 


No peaks 


5 C 


Bump-B, C 


C, B 


B 


No peaks 



rates at the end of region 3. Its fractional rms is higher in the hardest 
band (13-38 keV); 

(ii) a second peak is detected in region 2 (in three observations, 
see table[3) with properties markedly different from the first one. It 
has initially (in region 2) a frequency between 2.44 Hz and 2.84 Hz 
and a quality factor Q < 2. In regions 3 and 4 it shows a much less 
variable frequency than the first peak (slowly increasing from ~4 to 
7 Hz) and a Q > 2. Its fractional rms is higher in the hardest band 
(13-38 keV), is lower than that of the first peak and is independent 
on the centroid frequency (see Fig.|5]l- 

All five analyzed observations follow the average behaviour 
described above, with three main exceptions: 

(a) region 1, which correspond to the long hard intervals in the light 
curve of class /3 (see Fig.[T] upper panel), are not present in the two 
class fi observations (#1 and #4). Since this is the main macroscopic 
difference between the two classes, and since a detailed analysis of 
the long hard intervals of class /3 light curves is beyond the scope 
of this work, in the following of the paper we will consider the two 
classes as a single class. 

(b) No QPO is detected in the PDS from region 4 of observations #2 
(class /3) and #3 (class fj.). In region 4 of observation #3, only in the 
highest energy band, we find a non-significant QPO (a=2.4) with 



the following properties: v = 5.85±0.26 Hz, Q = 3.21±2.43, rms 
= (6.26±1.35)%. Both the rms amplitude value and the centroid 
frequency are consistent with those of the peak found in region 4 in 
the 36-103 keV energy band for observations #1, 4 and 5, but the 
uncertainty on the quality factor Q and the low significance (<3 a) 
does not allow us to claim the real detection of a QPO. 
(c) In region 2, we detected two significant peaks in observations 
# 2, 3 and 5, while just one significant peak in observations #2 
and 4. In region 2 of observation #4, we find a non-significant 
QPO (a=2.44) with the following properties (values from 14-35 
keV energy band): i/ = 3.46±0.18 Hz, Q = 2.44±1.03, rms = 
(3.68±0.93)%. Its properties are consistent with those of the low- 
frequency peak found in region 2 for observations #1,3 and 5, but 
the significance (<3 a) does not allow us to claim the real detection 
of a QPO. 

We will refer to the peak with centroid frequency ranging from 
~2 up to 15 Hz as type-C QPO and to the peak with 2.44-7 Hz 
frequency detected in regions 2, 3 and 4 as type-B QPO. This iden- 
tification will be motivated below. In Figure [5] we show the frac- 
tional rms amplitude of all the detected Bumps/QPOs as a function 
of their centroid frequency, for each of the three energy bands ana- 
lyzed. In a few cases QPOs were detected only in one or two energy 
ranges. In these cases in order to obtain a more reliable estimate of 
the peak amplitude also in the energy range where it was not signif- 
icant we fixed its centroid and FWHM to the values obtained from 
the energy range in which it was more significant. 

The points with large errors in frequency correspond to detec- 
tions where we need two Lorentzian components to obtain a satis- 
factory fit. For these points we plot the average of the two centroid 
frequencies with an error estimated as the sum between the two er- 
rors (we considered this the most conservative approach possible). 
In order to check whether these double peaks are due to the fre- 
quency variability of the QPO we tried to select smaller regions in 
the HHID. However we could not detect the QPO in such selec- 
tions, due to the low statistics of the resulting PDS. 

From Figure[3]we can relate the presence of the QPO-types to 
different zones of the source path in the HHID but we do not have 
any detailed information about the time evolution of the source po- 
sition in that diagram. We know that p and (i classes are charac- 
terized by fast and strong variability and frequent spectral transi- 
tions but we do not know precisely between which spectral states 
the transitions happen and when. A light curve where parts corre- 
sponding to different regions in the HHID are marked with different 
symbols can help us to clarify this point (see Figure [SJ. We chose 
a 300-second interval representative of the source behaviour. With 
the identifications Downward Triangles = Reg. 7, Circles = Reg. 
2, Stars = Reg. 3, Upward Triangles = Reg. 4 and Squares = Reg. 
5 we have the following regions sequence: 2345 345 532. Class /3 
observations are also characterized by the presence of long hard in- 
tervals (region 1), occasionally occurring when the source maps the 
hard tail in the bottom-right part of both the diagrams in Figure [2] 
always between region-2 intervals. If, using the state classification 
presented in Belloni et al. (2000), we make the identifications Reg. 
2 = C-state, Reg. 4 = A-state, Reg. 5 = B-state and we consider Reg. 
i as a Transition (T) state, we have the following state sequences: 
C T AB T AB T C. In this framework region 1 is the hardest part 
of the C-state. 



3.1 Phase lags 

Phase lags for all the detected peaks are reported in Tables |A6l and 
IA7l in appendix IaI In Figure [7] we show phase lags vs. frequency 
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Figure 5. Fractional rms of the detected Bumps/QPOs in all observations as function of their frequency for energy bands 2-5 keV (PCA channels 0-13, left 
panel), 5-13 keV (PCA channels 14-35, center panel) and 13-38 keV (PCA channels 36-103, right panel). Error bars are estimated at 1 ct confidence level. 
Arrows indicate 3 cr upper limits for type-C QPOs only. Filled symbols indicate type-C QPOs, open symbols type-B. 
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Figure 6. 300-second light curve (2-13 keV energy band) for observation # 5. Bin size is 0.125 sec. Every symbol is averaged on 2 seconds. Different symbols 
correspond to different regions in the HHID in Figure[3] 



for the type-C QPO in two observations. If the QPO was fitted with 
two Lorentzians we estimated the phase lag integrating on a peak 
centered on the centroid frequency of the Lorentzian with the high- 
est rms amphtude and with the FWHM of the same component. For 
observation # 1, a negative trend for increasing QPO frequencies is 
visible, although one point is completely out of this trend. For ob- 
servation # 5 the phase lag of the QPOs are always negative and 
consistent with being constant. 



In Figure [8] we show two examples of power density spec- 
tra (upper panels) and phase-lag spectra (bottom panels) from re- 
gions where a type-B QPO was detected. In the left panel (ob- 
servation #4), even thought the errors on phase lags are large, a 
turn towards positive lag values in correspondence of the QPO fre- 
quency is apparent. However, extracting the phase lags over the 
whole peak width (i.e. over the frequency range i/p ± FWHM/2) we 
obtain a negative value of the lags (although consistent with zero, 
= -0.07±0.05 rad). In the right panel (observation #1) we have 
a much better statistics and the phase lags are clearly negative ((f} 
= -0.11±0.007 rad) in correspondence of the type-B QPO (when 
integrated over the whole peak width). 



Filled symbols: obs, 1 
Empty symbols: obs 5 

Circles: C — type (regions 1 and 2) 
Squares: C — type (region 3) 



Frequency (Hz) 

Figure 7. Phase lags of type-C QPOs detected in two representative ob- 
servations. For each QPO we extracted the phase lag in a range cen- 
tered at the QPO peak frequency and corresponding to the width itself 
(Up ± FW H M /2). Errors bar on the X-axis are not shown for clarity. 
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Figure 8. Power density spectra and phase-lag spectra of two different 
observations/regions where a type-B QPO was detected. The dashed hne 
marks the QPO frequency, the dotted Hnes mark u it HWHM. Left panek 
type-B QPO detected averaging power spectra (121 points) in region 4 of 
observation # 5. Right panel: type-B QPO detected averaging power spectra 
(842 points) in region 3 of observation # 1 . 




Figure 9. Phase lags of the detected QPOs in all observations. For each 
QPO we extracted the phase lag in a range centered at the QPO peak fre- 
quency and coiTesponding to the width itself (Up ± FW H M/2). Errors 
bar on the X-axis are not shown for clarity. 



4 DISCUSSION AND CONCLUSIONS 

We analyzed 5 RXTE/PCA observations collected during the first 
two years of the mission. In the power density spectra of all five 
observations we detected several peaks which we identify with two 
different types of QPO already seen in many other BHCs: the type- 
C and type-B QPOs. This is the first identification of a type-B QPO 
in GRS 1915-1-105. To detect it, we looked in detail at spectral tran- 
sitions in observations characterized by a fast and intense variabil- 
ity. Spectral transitions in GRS 1915-1-105 are usually very fast, of- 
ten occurring on timescales of ~ seconds. This is at variance with 
most of other black-hole binaries in which spectral transitions are 
observed to last hours or days. In order to study the spectral transi- 
tions in GRS 1915-1-105 we performed an energy-dependent timing 
analysis by averaging power spectra on the pattern the source re- 
cursively tracks in the 3-dimensional hardness-hardness-intensity 
diagram. 

Applying this method, we found a type-B QPO in all five ob- 
servations. In all of them, we detected the type-B QPO together 
with the type-C, in region # 3 of the HHID (see Tab. [3]l- In three 
of them, we also detected a type-B QPO alone, in region # 4 of the 
HHID. Two of these observations belong to class fi and one to class 
/3, which excludes any relation between the class of variability and 
the presence of the type-B QPO. This means that the presence of 
the long hard intervals (which differentiate the /9-class from the /i- 
class light curves) does not influence the fast timing properties of 
the source outside these intervals. We could not find any property 
(as e.g. hardness, rate) correlated with the presence of the type-B 
QPO alone in region 4 of the HHID. In three observations we also 
found a type-B bump in region 2. No correlations were found nei- 
ther between the presence of the type-B bump and the type-B QPO 
in region 4 nor with the hardness and the count rate. 



behaviour very similar to that of type-C QPOs: its frequency is cor- 
related with the count rate and is inversely correlated with the hard- 
ness. At a certain hardness, this QPO disappears. A second type of 
QPO is also detected: as in the case of type-B QPOs, this second 
QPO appears in a narrow frequency range (often around r--^6 Hz) 
and in a limited range in hardness. 

Its frequency and quality factor appear to be slightly corre- 
lated with the count rate, particularly when at its lowest frequencies 
(2.44-2.84 Hz). We interpret this as the first evidence of an increase 
of the coherence of the type-B QPO from Q < 2 when at low fre- 
quencies to Q > 2 when reaching frequencies around ~6 Hz, pos- 
sibly suggesting the presence of a resonance at this frequency (see 
Casella et al. 2004). 

The combined evolution of the QPOs in GRS 1915-1-105 is 
strongly reminiscent of the known behaviour of type-C and type-B 
QPOs in BHCs (see Casella et al. 2005 and references therein). 

To verify this identification, we plot in Figure [5] the nns frac- 
tional variability of the detected QPOs as a function of their fre- 
quency for three energy bands. The two QPO types have a some- 
what similar energy dependency (being stronger at high energy) but 
they clearly show different behaviours in these diagrams. In each of 
the three panels, two well-identified groups of points are evident. A 
comparison of these two groups with Figure #3 of Casella et al. 
2004 helps to classify the observed QPO in GRS 1915-1-105: the 
first group is diagonally spread across the plots, covering the whole 
frequency range between ^ 2 and ~ 15 Hz and a large range in rms 
(particularly at high energies, see the right panel of Figure|5}. The 
second group is clustered both in frequency (between ~ 2.5 and 
~ 7 Hz) and in fractional rms. The observed behaviour is clearly 
consistent with that known to be typical of type-C and type-B QPOs 
in BHCs (see Casella et al. 2004, 2005). 



4.1 QPOs identification 

In order to identify the two types of QPO that we found in our data 
sets, we compare them with the known LFQPOs in BHCs. In partic- 
ular we first analyze their position and behaviour in the Hardness- 
Intensity diagram (right panel of Figure|2}. When the source moves 
through regions #1 and #2 up to region #3 the first QPO shows a 



4.2 Phase lags 

The association between QPO-type and phase lag in literature is 
not conclusive: although an average behaviour can be identified 
(see Tab. [T] Casella et al. 2005 and reference therein) there are a 
number of exceptions (see e.g. Belloni et al. 2005 and Homan et al. 
2005a). Nevertheless we performed a phase-lag analysis in order to 
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have a comprehensive view of the behaviour of the QPOs in GRS 
1915+105. 

On the basis of the analyzed data it is not possible to char- 
acterize unambiguously the phase lag behaviour of any of the two 
types of QPO we observe: lags appear to vary between different 
observations, although both types show in average negative values 
of phase lag (see Fig.|9]l. 

This can possibly due to the fast variability of the analyzed 
light curves: we extracted power spectra over time intervals 2 sec- 
onds long, without applying any procedure to detrend the variabil- 
ity on longer time scales. However, this variability is very strong 
(see Fig. [TJ, which results in a leakage at higher frequencies as 
strong as to actually dominate the phase-lag continuum. 



4.3 GRS 1915+105 as a "normal" source 

The type-B QPO was detected in a few BHCs (see Casella et al. 
2005), and associated to spectral transitions from HIMS to SIMS 
(Homan & Belloni 2005, Belloni et al. 2005 and reference therein). 
Some authors (Gallo et al 2003; Fender, Belloni & Gallo 2004) sug- 
gested a relation between these spectral transitions and the emis- 
sion and collimation of transient relativistic radio jets: if we con- 
sider the type-B QPO as the signature of radio jet emission, we 
have this "signature" also in the prototypical galactic jet source. 
As already pointed out by Belloni et al. (2005), the non-detection 
of a type-B QPO in GRS 1915+105 was rather interesting, espe- 
cially if you interpret the X-ray /radio correlation in this source and 
in other transients in the framework of the same model (Fender, 
Belloni & Gallo 2004). In GRS 1915+105 the association between 
X-ray and radio activity is well known (Pooley & Fender 1997, 
Mirabel et al. 1998, Fender & Belloni 2004 and references therein). 
Belloni et al. 2005 suggested that the elusiveness of this QPO (in 
GRS 1915+105) could be due to high velocity of the movement of 
the source through the HID. The result presented in this work thus 
strengthens the interpretation of GRS 1915+105 in the framework 
of the same model of other BHCs: a type-B QPO appears in cor- 
respondence of spectral transition from the B-state to the A-state 
(Transition, Region 3, Stars in Figure|6ll and when the source is in 
the A-state (Region 4, Upward Triangles in Figure [6]l. In the light 
curve in Figure|6]we see state oscillations CTAB TAB TAB that we 
can identify as fast passages between the HIMS and the SIMS of 
the other BHCs (Casella et al. 2004; Fender, Belloni & Gallo 2004). 
The type-B QPO appears also in correspondence of the transition 
from the B-state to the C-state (Squares-Stars-Circles in Figure|6ll, 
therefore not in an oscillation event but in a transition from a soft 
to a hard state. 

Furthermore, GRS 1915+105 is at present the heaviest black 
hole (although uncertainties on BH mass values are rather large, see 
McClintock & Remillard 2006 and references therein) in which a 
type-B QPO has been observed at 6 Hz (when its quality factor 
Q > 2). This strengthens the already known independence of this 
type of QPO on the mass of the compact object (see Casella et al. 
2005). 

4.4 X-ray/radio association 

Klein-Wolt et al. (2002) made a detailed study of simultaneous 
radio/X-ray observations of GRS 1915+105, focusing in particular 
on radio oscillation events, and found that {3 class observations are 
usually associated to strong radio oscillations, while class obser- 
vations are usually associated to a weak steady radio activity. These 



authors point out that the long hard intervals (which are the main 
macroscopic difference between the two variability classes) appear 
to be directly linked to the production of radio oscillations. We 
found the type-B QPO in both classes 13 and ^l. We could not find 
any timing property (on time scales shorter than 2 seconds) clearly 
differentiating the two classes. This apparently weakens the link 
between the presence of the type-B QPO and radio activity. Un- 
fortunately, the lack of radio coverage during the analyzed RXTE 
observations does not allow any conclusive result. In order to clar- 
ify this issue it will be fundamental to extend our analysis to RXTE 
observations for which a radio coverage is available. 
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APPENDIX A: FIT RESULTS AND PHASE LAGS 

Here we report five tables containing the fit results for the anal- 
ysed observations and two tables with the results of the phase lag 
analysis. 



Table Al. Fit results for Obs. #1. We accounted for the Poissonian noise 
adding an additive constant to the model (see f|2). In the first column the 
bold number indicates the region, the second number indicates the selec- 
tion within the region and it increases following the arrows in figure |2] 
For peaks fitted with two Lorentzians, both the components are reported. 
If the difference between the negative and the positive error was smaller 
than 0.1, we reported just the value of the largest error (in absolute value). 
With "Fixed fit" we refer to those cases when QPOs were detected only in 
one or two energy ranges: in these cases we fixed the centroid frequency 
and the FWHM to the values obtained from the energy range in which 
the peak was more significant (see [[3). 

Observation #1 



Sel. 


Type 


Freq. (Hz) 


FWHM 
2-5 keV 


Sign, (cr) 


rms (%) 


2-1 


B 


2.75±0.16 


, 79+0.42 
'■'^-0.36 


4.68 


4.18±0.45 


2-1 


C 


9.35±0.33 


'•°'-0.68 


3.10 


2.04±0.39 


3-1 


B 


3.99±0.07 


2.21±0.21 


13.26 


3.46±0.13 


3-1 
3-2 


C 
B 


10.92±0.32 
5.04±0.06 


43^+0.88 
2.39±0.21 


7.00 
16.13 


2.05±0.16 
3.15±0.10 


3-2 


C 


1 ^ nQ+0.48 




3.45 


1.25±0.23 


3-3 


B 


5.65±0.09 


2.70±0.27 


14.00 


3.33±0.12 


3-3 


C 


12.43 


4.39 


Fixed fit 


0.91±0.20 


3-4 


B 


6.01±0.16 


2.68±0.49 


8.53 


2.06±0.13 


3-4 


C 


14.53 


2.66 


Upp. lim.'^ 


0.70 


3-5 


B 


5.94 


5.51 


Upp. lim.'^ 


1.47 


4-1 


B 


6.39±0.20 


1 99+0-91 

^-'^-0.62 


4.48 


3.07±0.41 


5-13 keV 


2-1 


B 


2.75 


1.72 


Fixed fit 


5.23±0.52 


2-1 


C 


8.67±0.27 


5 42+0.95 
-'-^^-0.83 


7.87 


7.02±0.45 


3-1 


B 


4.27±0.04 


1.84±0.10 


23.23 


5.60±0.12 


3-1 


C 


10.67±0.10 


3.95±0.27 


19.09 


3.70±0.10 


3-2 


B-lst 


4.81d=0.08 


1.39±0.18 


5.56 


a 


3-2 


B-2nd 


5.72±0.06 


1.03±0.11 


5.01 


5.23+0.50*= 


3-2 
3-3 


C 
B-lst 


12.80±0.40 
5.28±0.12 


o 92+1-56 

1.81±0.20 


6.74 
8.08 


2.40+0.18 

_a 


3-3 


B-2nd 


6.34±0.06 


0.82±0.14 


4.48 


6.38+0.52*= 


3-3 
3-4 


C 
B 


12.43±0.58 
6.10±0.05 


4 ,g+2.26 

^■-'^-1.52 

2.32±0.15 


3.56 
21.58 


1.58+0.27 
4.20+0.10 


3-4 


C 


14.53±0.35 


T <;/;+l-10 
2.66I0 72 


4.27 


1.16+0.16 


3-5 


B 


. Q4+0.37 
^■"-0.50 


5 51+1-12 


6.39 


3.42+0.35 


4-1 


B 


6.78±0.08 


(,(-,+0.80 


5.43 


T ^^ + 1.55 
3-66I0.34 


13-38 keV 


2-1 


B 


2.74 


1.72 


Upp. lim.'* 


4.38 


2- 1 

3- 1 
3-1 


C 
B 
C 


9.08±0.30 
4.27 
11.24±0.10 


4 01+0-96 
^■"1-0.84 

1.84 
3.53±0.39 


7.43 
Fixed fit 
10.91 


17.25+1.20 

4 92 + 1-35 
^-^^-0.50 

12.32+0.56 


3-2 
3-2 
3-2 
3-3 


B 

C 
C<= 
B 


5.77±0.15 
12.65±0.13 

17.88+277 

5^65 


, 07+0.54 

'•-"-0.42 
7 77+0.54 
■^■"-0.48 

19.96t?:°^ 
2.70 


3.93 
5.19 
3.76 
Fixed fit 


4.80+0.69 
8.70+0.87 

6.82+0^50 


3-3 


C 


14.53 


2.66 


Fixed fit 


2.89+0.53 


4-1 


B 


6.78 


0.63 


Fixed fit 


7.83+0.67 



^ Peak fitted with two Lorentzians. The value of the rms is reported at the 
next line 

*= Sum of the tins of the two Lorentzians fitting this peak 
Harmonic peak 

Estimated at Scr confidence level 
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Table A3. Fit results for Obs. #3. In the first column the bold number indicates the region, the second number indicates the selection within the region and it 
increases following the arrows in figure|2] For peaks fitted with two Lorentzians. both the components are reported. If the difference between the negative and 
the positive eiTor was smaller than 0.1, we reported just the value of the largest error (in absolute value). With "Fixed fit" we refer to those cases when QPOs 
were detected only in one or two energy ranges: in these cases we fixed the centroid frequency and the FWHM to the values obtained from the energy range 
in which the peak was more significant (see ij3). 



Observation #3 


Sel. 


Type 


Freq. (Hz) 


FWHM 
2-5 keV 


Sign, (cr) 


rms (%) 


1-1 


C 


2.72 


2.29 


Fixed fit 


10.12±0.35 


1-2 


C 


3.58±0.09 


2.41±0.28 


6.57 


6.42±0.53 


1-3 


c 


5.80±0.10 


1.93±0.28 


8.66 


4.36±0.26 


1-3 




13 41+0.71 


5 6n+'''2* 

■'■°"-1.98 


3.74 


240+0.45 

^■™-0.33 


2-1 


B 


2.44±0.28 


2.24±0.49 


3.06 


2 95+0.68 


2-1 


C 


7.57±0.07 


2.28±0.24 


21.08 


3.14±0.13 


2- 1 

3- 1 


B 


15 91+3.26 
5.03±0.33 


19 98+0.02 


5.01 
4.00 


2 01 +"•3° 
■^•"^ -0.20 
1 77+0.33 
^- "-0.22 


3-1 


C 


1I.00±0.21 




6.94 


1.84±0.13 


3-2 


B 


5.91±0.09 


0.71±0.18 


5.46 


1.47±0.13 


3-2 


C 


12.59±0.27 


3.69±0.84 


6.79 


1.70±0.15 


3-3 


B 


5.53 


0.75 


Fixed fit 


0.41±0.28 


3-4=- 


B 


5.37 


1.45 


Fixed fit 


0.51±0.10 


5-13 keV 


1-1 


C 


2.72±0.14 


4.01±0.23 


15.77 


18.04±0.58 


1-2 


C 


3.35±0.08 


2.80±0.22 


8.23 


14 13+103 


1-3 


C 


5.71±0.09 


2.13±0.21 


9.09 


8.61±0.47 


1-3 


C<: 


13.90±0.54 


10.36ti;23 


9.81 


5.58dz0.31 


2-1 


B 


2.64±0.21 


2.24±0.38 


3.50 


4 4Q+0.76 
^■^"-0.63 


2-1 


C 


7.65±0.04 


1.94±0.10 


24.82 


5.74±0.12 


2-1 




14.S5_o 92 


1 5 00+1-82 


9.05 


3.88±0.25 


3-1 


B-lst 


5.77±0.03 


<0.5f 


6.93 


_b 


3-1 


B-2nd 


4.17±0.12 


2.31±0.41 


5.42 


4.04±0.43'l 


3-1 


C 


11.17±0.12 


3.89±0.29 


15.98 


3.13±0.10 


3-1 


_e 




Q "37+3.35 
'■•"-2.23 


3.80 


1.27±0.20 


3-2 


B-lst 


5.17±o!o6 


1.00±0.14 


13.78 


_b 


3-2 


B-2nd 


6.25±0.02 


<0.5f 


8.23 


3.97±0.18'* 


3-2 


C 


13.19±0.16 


3.66±0.45 


13.39 


2.41±0.11 


3-3 


B 


5.53±0.05 


0.75±0.23 


8.14 


1.99±0.13 


3-4^ 


B 


5.36±0.16 


1 45+0.45 


3.90 


0.78±0.10 


13-38 keV 


1-1 


C 


2.80±0.22 


3.42±0.36 


8.13 




1-2 


C 


3.20±0.10 


3.18±0.20 


11.42 


21.77±0.96 


1-3 


C 


6.15±0.07 


1.29±0.30 


5.95 


10.61 ±0.89 


2-1 


B 


2.64 


2.24 


Fixed fit 


5.29±0.69 


2-1 


C 


7.71±0.06 


1.58±0.14 


9.84 


11.65±0.59 


2- 1 

3- 1 


B 


17.72t-^? 
5.29±0.21 




3.55 
4.32 


7 17+0.92 
'■•"-1.05 

4 65+0.75 

^■°^-0.54 


3-1 
3-1 


C 


11.35±0.11 
20.371^;^° 


4.21±0.44 
,9 95+0.05 


7.80 
3.93 


ii.o5ir7? 


3-2 


B 


5.38±0!08 


1.47±0.19 


10.34 


8.39±0.43 


3-2 


C 


12.57±0.19 


2.98±0.54 


4.02 


7 70 + 0.86 
'■'•'-0.96 


3-3 
3-4^ 


B 
B 


4.72±0.47 

401+0.33 
^■^^-0.53 


T. Q7+l.'''8 
-'■'^'-1.07 


4.14 
3.44 


5 47+0.88 
•^■^'-0.67 
9 ,7+0.70 
-^■•"-0.35 



Poissonian noise fitted adding a constant to the model, since the Poissonian noise subtraction left flat residuals at high frequencies (see ijJJ 
^ Peak fitted with two Lorentzians. The value of the rms is reported at the next line 
Harmonic peak 

Sum of the rms of the two Lorentzians fitting this peak 
" High frequency QPO, not considered in our analysis (see Belloni et al. 2006) 
^ Our resolution is 0.5 s (we extracted power density spectra on 2-second intervals, see |j2j 
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Table A2. Fit results for Obs. #2. See caption of Table |A3l for an explana- 
tion. 



Table A4. Fit results for Obs. #4. See caption of Table |A3l for an explana- 
tion. 



Observation #2 



Observation #4 



Sel. 


Type 


Freq. (Hz) 


FWHM 
2-5 keV 


Sign, (cr) 


rms (%) 


Sel. 


Type 


Freq. (Hz) 


2-5 keV 


Sign, (a) 


rms (%) 


J.- i 




i.oomu.Lfo 


1 9«;_i_n 1 Q 
1 .zjmu. Ly 


0.40 


'^ OA -I-H Af\ 




\^ 


1 ^7 
iU. JZ 


J. "40 


r'lxea lit 


7 1 q+0 77 

z. 1 ymu.z / 


l-i 




A QQ_l_n A 1 


J.41itU.Oo 
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C 
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11 Qi +n 7 1 
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•'■'-'-0.55 




1 c 1 +n 1 ^ 


1-3 


c 


6.06±0.13 


1.74±0.33 
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B 
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Q 


7 72it0 13 




7 06 
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\^ 
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Q 
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7.1 1 
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T) 

13 
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B 
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Fixed fit 
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R 
13 


O.O^ZtU. i J 


7 nf»+n 47 

Z.UDlllU.4Z 


J.Z4 


^ 1 9-1-0 n9 

J. 1 / itu.uz 


3-2 


C 


12.72 


5.51 


Upp. lim.<^ 


4-1 


B'^ 






3.70 


7 -,0+0.46 

' _0 34 








5-13 keV 












5-13 keV 






1-1 


C 


2.96±0.58 


1.43±0.13 


11.34 


13.36±0.60 


2-1 


C 


9.92±0.21 


3.09±0.50 


6.70 


4.80±0.36 


1-1 




5.25±0.61 


5.44±0.61 


5.16 


o 7, + 1.10 
8-2i__o.gi 


3-1 


B 


4.27±0.04 


2.01±0.10 


24.77 


6.20±0.13 


1-2 


C 


4.18±0.06 


2.15±0.21 


11.76 


9.73±0.47 


3-1 


C 


11.45±0.11 


4.34±0.39 


18.25 


3.71±0.12 


1-3 


C 


6.06±0.07 


1.25±0.19 


9.75 


5.31±0.31 


3-2 


B-lst 


4.82±0.06 


1.38±0.12 


10.15 


_a 


2-1 


C 


7.91±0.06 


1.74±0.16 


25.56 


5 29+0-20 
-^-^'-0.10 

T ,t; + 0.58 

•'--^-'-0.39 

3.62±0.13 


3-2 


B-2nd 


5.79±0.05 


0.66±0.12 


4.75 


5.93±0.39'' 


2-1 
2-2 


c 


1-'-^ -2.18 

9.81±0.08 


15.48t^-5i 
2.41±0.26 


4.68 
14.20 


3-2 
3-3 


C 
B 


12.52±0.57 
5.79±0.05 


q q4.+l-85 
^•^+-1.56 

1.61±0.10 


7.49 
22.08 


2.54±0.19 
5.87±0.I4 


2-2 




21.15±0.92 


7 ,,+3.80 
'-Jl_2.32 
<0.5<= 


4.68 


, <:q+0.32 
1-08_Q ig 

_a 


3-3 


C 


13.63 


6.72 


Fixed fit 


0.97±0.27 


3-1 


B-lst 


5.75±0.01 


8.07 


3-4 


B 


5.61±0.05 


1.74±0.16 


16.41 


4.12±0.13 


3-1 


B-2nd 


4.57±0.10 


1.39±0.23 


7.09 


4.59ib0.33b 


3-4 


B<= 


11 n+o-62 




3.89 


1.63±0.27 


3-1 


C 


12.78±0.12 


2.1±50.34 


9.09 


2.23ib0.13 


4-1 


B 


6.79±0.07 


1.66±0.22 


10.79 


6.48±0.30 


3-2 


B 


5.10±0.10 


0.98±0.21 

, 77+0.48 
^■'^-0.38 


7.08 


3.80ib0.27 


4-1 


B<^ 






6.68 


5.54±0.34 


3-2 


C 


11.75±0.18 


5.47 


1.83ib0.19 








13-38 keV 







13-38 keV 



1-1 


C 


2.89±0.06 


1.67±0.21 


6.56 


20.17±1.57 


1-2 


C 


4.01±0.17 


3.54±0.43 


6.06 


17.62±1.46 


1-3 


c 


6.55±0.10 


1.59±0.35 


6.02 


10.24±0.85 


1-3 




IT 1Q+2.26 
1^-1^-2.67 


1 1 ,'1+3.66 
ll.JJ_3 26 


3.37 


7.85il:fg 


2-1 


c 


7.97±0.13 


1.72±0.39 


4.10 


10.74ll:« 


2-2 


c 


9.84±0.14 


1.73±0.35 


5.19 


8.47lg:|^ 


3-1 


B 


5.16±0.26 


l-^'-0.55 


5.45 


7.04±0.71 


3-1 


C 


12.94±0.18 


5.05±0.33 


10.36 


9.92±0.48 


3-2 


B 


4.74±0.245 


1 36+0-56 
l-^°-0.21 


6.07 


o ,0+1-40 

**-l'J-o.67 


3-2 


C 


,7 77+0.43 
lz.;z_Q 29 


, j-,+0.86 
-'-'1-1.06 


7.07 


n Q7 + 0.52 
^-^^-0.71 



2-1 


C 


10.32±0.19 


3 46+"-''^ 
^-™-0.53 


9.84 


13.68±0.71 


3-1 


B 


4.27 


2.01 


Fixed fit 


4.02±0.99 


3-1 
3-2 
3-2 


C 
B 
C 


11.67±0.13 
4.62±0.3I 
I2.67±0.19 


2.82±0.44 
7 ,,+0.96 

Z.JJ_Q yg 

5.91±0.76 


6.90 
3.46 
13.31 


9.49±0.69 

5- 70 + 1.00 

-'■^°-0.77 

11.16±0.45 


3-3 
3-3 

3- 4 

4- 1 


B 
C 
B 
B 


5.85±0.30 
13.63±0.54 

^41+0-62 
-'-^1-0.79 
6.83±0.14 


2.39±0.63 

g 77+2.66 
°-'^-1.67 
gQ,+3.50 
D-Ul-2.24 
, ^7+0. 45 
l-'^-0.35 


4.01 
6.39 
4.66 
5.28 


7 7A+0.65 
'■^l'-0.90 
0-7+1.68 
<*-^'-0.94 

5 74+1-54 
-"■'^-0.65 
,0 76+1-78 



^ Peak fitted with two Lorentzians. The value of the rms is reported at the 
next line 

^ Sum of the rms of the two Lorentzians fitting this peak 
Harmonic peak 

Estimated at 3a" confidence level 
° Our resolution is 0.5 s (we extracted power density spectra on 2-second 
intervals, see ^ 



^ Peak fitted with two Lorentzians. The value of the rms is reported at the 
next line 

'' Sum of the rms of the two Lorentzians fitting this peak 
Harmonic peak 
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Table AS. Fit results for Obs. #5. See caption of Table |A3l for an explana- 
tion. 

Observation #5 



Sel. 


Type 


Freq. (Hz) 


FWHM 

A- J Kc V 


Sign. (<t) 


nns (%) 


1-1 


C 


2.53±0.22 


3.92±0.41 


7.99 


^■^"-0.56 


1-2 


C 


3.17±0.20 


4.61±0.34 


9.50 


8.44±0.45 


1-3 


c 


6.42±0.05 


1.26±0.13 


14.33 


3.28±0.13 


1-3 




14.03±0.65 


O cn + 2.38 

o-Jo_1.74 


6.46 


2.50±0.27 


2-1 


B 


2.75±0.10 


2.50±0.39 


13.13 


2.87±0.13 


2-1 


C 


8.40±0.16 


3.14±0.40 


4.42 


2.53±0.38 


3-1 


B 




5 ,4+0.83 
-^■^^-0.73 


4.80 


3.41±0.36 


3-1 


C 


11.88±0.15 


, 7, +0.66 
^■'^-0.50 
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3-2 


B 
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2.59±0.39 
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3.09±0.23 


3-2 


C 


12.85±0.33 
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5.66 
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3-3 


B 
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3.28 


2.26±0.41 


3-3 


C 


15.31 


6.22 


Fixed fit 
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B 
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5-13 keV 
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1-2 


C 
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3.05±0.44 
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1-3 
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5.87±0.24 


2-1 


B 
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2-1 


C 
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2-1 




18.43±0.59 
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3-1 
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4.40±0.12 


2.85±0.15 


10.30 


_a 


3-1 
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6.25±0.005 


<0.5° 


7.94 


5.54±0.29'' 


3-1 


C 


12.08±0.07 


2.27±0.22 


15.28 


2.49±0.08 


3-2 
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5.07±0.08 
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6.57 


_a 
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6.42±0.10 
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C 
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2.28±0.19 
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B 
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1.85±0.55 
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13-38 keV 


1-1 


C 


2.80±0.11 
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1-2 


C 
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3.84±0.35 


5.81 
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6.69±0.04 


1.76±0.14 
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2.84 
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1.83 


2-1 
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10.17±0.17 


2.67±0.33 


5.81 


_a 


2-1 


C-2nd 


8.34±0.09 


1.40±0.29 


4.31 


11.92+i-?^b 


3-1 


B 


5.14±0.29 


2 ^4+0.68 


4.43 


"■^^-0.71 


3-1 


C 


11. 81 ±0.09 


2.75±0.34 


8.80 


9.18±0.55 


3-2 


B 


6.04±0.11 


0.97±0.26 


4.57 


6.58±0.72 


3-2 


C 


13.56±0.23 


4 27+0.71 
^■^'-0.61 


10.04 


9.42±0.49 


3-3 


B 


5.35±0.18 


9 99+0.63 
^■^^-0.50 


5.53 


7Q6+0.70 

'■"°-0.64 


3-3 


C 


17 (S5+-'^'^^ 
^ '■°-'-i.64 


1 1 -7/S+4.12 


4.86 


'■-'"-0.77 


4-1 


B 


5.96 


0.76 


Fixed fit 


4 02+0-87 
^■"^-1.01 



Table A6. Phase lags for each component detected in the fits for obser- 
vations #1,2 and 3, integrated in a range centered on the peak centroid 
frequency with a width equal to the FWHM of the peak itself. Centroid 
frequency values were taken from fits in the 5-13 keV energy band. For 
peaks fitted with two Lorentzians, we calculated the lag integrating in a 
range centered on the component with the highest fractional nns (see ||2). 

Observation #1 

Sel. Type Freq. (Hz) Phase lag (rad) 



2-1 


B 


2.75 


-0.042±0.020 


2-1 


C 


8.67±0.27 


0.051±0.019 


3-1 


B 


4.27±0.04 


-0.110±0.007 


3-1 


C 


10.67±0.10 


0.040±0.007 


3-2 


B=» 


4.81±0.08 


-0.142±0.008 


3-2 


C 


12.80±0.40 


0.001±0.005 


3-3 


ga 


5.28±0.12 


-0.089±0.001 


3-3 


c 


12.43±0.58 


-0.019±0.001 


3-4 


B 


6.10±0.05 


-0.025±0.010 


3-4 


C 


14.53±0.35 


-0.008±0.014 


3-5 


B 


5 94+0.37 


-0.019±0.0I3 


4-1 


B 


6.78±0.08 


-0.095±0.031 


Observation #2 


1-1 


C 


2.96±0.58 


-0.018±0.022 


1-1 


Cd 


5.25±0.61 


-0.023±0.020 


1-2 


C 


4.18±0.06 


-0.120±0.020 


1-3 


C 


6.06±0.07 


-0.271±0.040 


2-1 


C 


7.9 1 ±0.06 


-0.293±0.055 


2-1 


Cd 


1 5 17+1.35 
13.1 /_2. 18 


-0.285±0.046 


2-2 


C 


9.81±0.08 


-0.115±0.044 


2-2 


C 


21.15±0.92 


-0.290±0.155 


3-1 




4.57±0.10 


-0.052±0.026 


3-1 


C 


12.78±0.12 


-0.123±0.066 


3-2 


B 


5.10±0.10 


0.140±0.060 


3-2 


C 


11.75±0.18 


-0.263±0.103 


Observation #3 


1-1 


C 


2.72±0.14 


-0.122±0.025 


1-2 


C 


3.35±0.08 


-0.093±0.013 


1-3 


C 


5.71±0.09 


-0.116±0.028 


1-3 


Cd 


13.90±0.54 


-0.031 ±0.053 


2-1 


B 


2.64±0.21 


-0.128±0.020 


2-1 


C 


7.65±0.04 


-0.189±0.024 


2-1 


Cd 




-0.202±0.021 


3-1 


ga 


4.17±0.12 


-0.101±0.019 


3-1 


C 


11.17±0.11 


-0.115±0.029 


3-1 


.b 


25.00tij!] 


-0.168±0.194 


3-2 


B 


5.17±0.06 


0.080±0.040 


3-2 


C 


I3.19±0.16 


-0.080±0.043 


3-3 


B 


5.53±0.05 


-0.282±0.120 


3-4<: 


B 


5.36±0.16 


0.008±0.060 



^ Peak fitted with two Lorentzians. The value of the nns is reported at the 
next line 

^ Sum of the rms of the two Lorentzians fitting this peak 
Harmonic peak 

Estimated at 3a confidence level 
° Our resolution is 0.5 s (we extracted power density spectra on 2-second 
intervals, see [|2) 



^ Peak fitted with two Lorentzians 
^ High frequency QPO 

Poissonian noise fitted adding a constant to the model 

Hai'monic peak 



A transient QPO from GRS 1915+105 

Table A7. Phase lags for each component detected in the fits for observa- 
tions # 4 and 5. See caption of Table|A6]for an explanation. 



Observation #4 



Sel. 


Type 


Freq. (Hz) 


Phase lag (rad) 


2-1 


C 


9.92±0.21 


0.021±0.102 


3-1 


B 


4.27±0.04 


-0.262±0.024 


3-1 


C 


11.45±0.11 


0.0±0.052 


3-2 


ga 


4.82±0.06 


-0.313±0.024 


3-2 


c 


12.52±0.57 


-0.188±0.046 


3-3 


B 


5.79±0.05 


-0.127±0.038 


3-3 


C 


13.63*^ 


-0.087±0.102 


3-4 


B 


5.61±0.05 


-0.122±0.041 


3- 4 

4- 1 


B'^ 
B 


11 13+062 
6.79±0.07 


-0.160±0.()96 
-0.074±0.048 


4-1 


B'^ 




-0.171 ±0.040 


Observation #5 


1-1 




2.69±0.27 


-0.130±0.023 


1-2 


C 


3.55±0.21 


-0.116±0.018 


1-3 


c 


6.46±0.03 


-0.1 25 ±0.024 


1- 3 

2- 1 


B 


12.86l»;22 
2.84±0.12 


-0.060±0.025 
-0.102±0.014 


2-1 


C 


8.70±0.05 


-0.061±0.019 


2-1 




18.43±0.59 


-0.057±0.064 


3-1 


ga 


4.40±0.12 


-0.140±0.015 


3-1 


c 


12.08±0.07 


0.033±0.046 


3-2 


ga 


5.07±0.08 


-0.267±0.036 


3-2 


c 


13.13±0.25 


-0.104±0.052 


3-3 


B 


5.40±0.07 


-0.290±0.050 


3-3 


C 


15.31±0.46 


-0.029±0.073 


4-1 


B 


e.lOitO.lO 


-0.0±0.044 



^ Peak fitted with two Lorentzians 

Fixed fit (see caption of table lATI for an explanation) 
Harmonic peak 



